Annealed disorder, rare regions, and local moments: A novel mechanism for 

metal-insulator transitions 
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It is shown that for noninteracting electron systems, annealed magnetic disorder leads to a new 
mechanism, and a new universality class, for a metal-insulator transition. The transition is driven 
by a vanishing of the thermodynamic density susceptibility rather than by localization effects. The 
critical behavior in d = 2 + e dimensions is determined, and the underlying physics is discussed. 
It is further argued that annealed magnetic disorder, in addition to underlying quenched disorder, 
describes local magnetic moments in electronic systems. 



PACS numbers: 71.30.+h; 64.60.Ak 

Metal-insulator transitions (MITs) remain a fascinat- 
ing and only incompletely miderstood phenomenon [y. 
Conceptually, one distinguishes between Anderson tran- 
sitions in models of noninteracting electrons, and Mott- 
Hubbard transitions of clean, interacting electrons. At 
the former, the electronic charge diffusivity D is driven 
to zero by quenched, or frozen-in, disorder, while the 
thermodynamic properties do not show critical behavior. 
At the latter, the thermodynamic density susceptibility 
dn/dfj, vanishes due to electron-electron interaction ef- 
fects. In either case, the conductivity a = {dn/d^)D 
vanishes at the MIT. In many real systems both quenched 
disorder and interactions are present, which makes a the- 
oretical understanding of the resulting Anderson-Mott 
transition very difficult. One particular complication 
is provided by the presence of magnetic local moments 
(LMs) in such systems. There is much experimental evi- 
dence for LMs [pi, and their formation has been studied 
theoretically |0] , but no existing theory can describe their 
interplay with the transport properties near the MIT pi . 
Another complication is the possible presence of annealed 
disorder, which is in thermal equilibrium with the rest of 
the system and hence involves disorder averaging of the 
partition function. This is in contrast to quenched disor- 
der that requires an averaging of the free energy, which 
is usually done by means of the replica trick Q . 

In this Letter we make two contributions to the MIT 
problem. (1) We show that annealed disorder leads to 
a MIT that belongs to none of the previously studied 
classes. It is driven by a vanishing dn/dpL and thus re- 
sembles a Mott-Hubbard transition, even if no correlation 
effects are explicitly considered. (2) We propose a mech- 
anism by which additional annealed disorder is generi- 
cally self-generated in quenched disordered systems, and 
we argue that a type of LMs can be described in terms of 
it. We further develop a method for incorporating these 
'annealed LMs' into a transport theory. 

Let us start by considering Wegner's nonlinear sigma- 
model (NLcrM) for noninteracting electrons with non- 



magnetic quenched disorder. The action reads 

^4=-^ /dxtr {\7Q{^)f + 2H /"rfxtr (OQ(x)). (1) 

Here Q(x) is a matrix field that comprises two fermionic 
degrees of freedom. Accordingly, Q carries two Matsub- 
ara frequency indices n and m, and two replica indices 
a and fi to deal with the quenched disorder. The ma- 
trix elements Q^m ^-^^ spin-quaternion valued to allow for 
particle-hole and spin degrees of freedom. It is convenient 
to expand them in a basis r^ ® Si (r, i = 0, 1, 2, 3) where 
To = So is the 2x2 unit matrix, and ti^2,3 = ^ si,2,3 = 
— icri,2,3, with Uj the Pauli matrices. For simplicity, we 
will ignore the particle-particle or Cooper channel, which 
amounts to dropping ti and T2 from the spin-quaternion 
basis Q|. Q is subject to the constraints <5^(x) = 1, and 
tr(5(x) = 0. 17"^j = dnmSap^n {tq (g) So) is a frequency 
matrix with 51„ = 27rTn a bosonic Matsubara frequency 
and T the temperature. G is a measure of the disorder 
that is proportional to the bare resistivity, and the fre- 
quency coupling H is proportional to the bare density 
of states at the Fermi level, tr denotes a trace over all 
discrete degrees of freedom that are not shown explicitly. 

The properties of this model are well known |||,^,y. 
The bare action describes diffusive electrons, with D = 
1/GH the diffusion coefficient. Under renormalization D 
decreases with increasing disorder until a MIT is reached 
at a critical disorder value. The critical behavior is known 
in an e-expansion about the lower critical dimension d = 
2. In the absence of the Cooper channel, the MIT appears 
only at two-loop order at a critical disorder strength of 
0{-\/e). H , which determines the specific heat coefficient, 
the spin susceptibility, and dn/dfj,, is uncritical, which 
makes this MIT an Anderson transition. 

Now we add magnetic annealed disorder to the model. 
Since our general results are independent of its origin, 
we ffi'st proceed without specifying it. Annealed disorder 
implies that the Q in the resulting terms all carry the 
same replica index H]; otherwise, the functional form of 



the resulting additional terms in the action can be taken 
from Ref. Q. We obtain two additional terms, viz. 

TAf r "^ r 



- tr (Q""(x))^ 
and Ay^(2) ^ Ay^(2,^) + Ay^'^.t)^ ^^ere 



(2a) 



A^(^.) = ^ E E E 

a^0 nm r=0,3 

x[tr(rt®s,)Q^"„(x)] 



dx[tr(r,®s,)g:^(x)] 



(2b) 



with j — s for i = (spin-singlet), and j = i for i = 1, 2, 3 
(spin-triplet). A^*^^^ arises from the need to absorb the 
scattering rate due to the annealed disorder in G. The 
coupling constants Mi, M2, and -M| are related to the 
strength of the magnetic annealed disorder. The factor 
of T appears naturally in front of any annealed disorder 
term, a crucial point that we will come back to later. 

The action ^+A^'^) + Ayl(^) can be analyzed by stan- 
dard means. Note that the mass terms in Eqs. (||) are 
proportional to temperature, making them quite differ- 
ent from conventional masses due to quenched disorder. 
In many respects, they are similar to electron-electron 
interaction terms in a Q-field theory formalism [Q. We 
denote the renormalized coupling constants that corre- 
spond to G, H, Ml, and A/^' by g, h, mi, and 7TI2' , and 
define 612 — mi'2/h. The renormalization group (RG) 
flow equations to one-loop order are 
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-e.g + .g2(5^+3(5* -3(5i) , 
:-/ig(<5^ + 3<5* - 3<5i) , 

-g[-A5l + 5i{5^2+i5\) + {5^2-S\?] 
-g[^5l + ^5i{6'2-'i5\)~i5'2{6l-5\)\ 



-i=g[i6l-5i{25i + 5^)-{5^f 



dl 
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(3a) 
(3b) 
(3c) 
(3d) 

(3e) 



where I — \nb with b the RG length scale factor. Besides 
unstable fixed points (FPs), there is a line of critical fixed 
points (FPs) {g*,h*,SlS'2*,Si*) = {e/AS*2,0,0, 6*2,6*2) 
that correspond to an MIT (all of these FPs belong to 
the same universality class). Linearization about any of 
these FPs yields one relevant eigenvalue A^ = e + O(e^) 
that determines the correlation length exponent j/ — 
1/Ag, one marginal eigenvalue that corresponds to mov- 
ing along the line of FPs, and two irrelevant eigenvalues 
equal to — e + O(e^). The anomalous dimension of h is 
K = — e + O(e^). In addition, the critical behavior of the 



single-particle density of states (DOS), N, at the Fermi 
level can be obtained from the wavefunction renormaliza- 
tion. Choosing the critical exponent of the DOS, /3, the 
correlation length exponent v, and the dynamical critical 
exponent z = d + n as independent exponents, we find 

u^l/e + 0(1) , (3^e + 0(6^) , z = 2 + 0(6^) . (4) 

For the conductivity exponent we find s — i/e = l + 0{e), 
and dn/d^, the spin susceptibility Xs, and the specific 
heat coefficient 7 = Cv jT , which we collectively denote 
by Xi all vanish with a critical exponent determined by 
K. The diffusion coefficient, on the other hand, has no 
anomalous dimension and thus is uncritical to one-loop 
order, as can be seen from Eqs. ( pa] , pq ). With t the 
dimensionless distance from the critical point at T = 0, 
and E the energy, we can summarize the critical behavior 
of these quantities by the homogeneity laws 



x{t,T) = b'^x{tb'/'',T¥) , 
N{t, T, E) ^ b-'^^^Nitb^^", Tb^Eb"") 
a{t,T)=b-'/''a{tb^/'',Tb'-) , 
D{t,T) = b-^''''+''^D{tb^''',TV) 



(5a) 
(5b) 
(5c) 
(5d) 



We conclude that the MIT is driven by the vanishing 
of dn/dfi, and therefore is qualitatively different from 
the localization transition that is found in the absence of 
annealed disorder. Indeed, putting Mi = M2 = M^ — 
we find that all thermodynamic anomalies disappear, as 
does the one-loop correction to g. At two-loop order, one 
finds instead a MIT of Anderson type ||l| . 

We now turn to a specific realization, via local mag- 
netic moments, of the annealed disorder that leads to the 
striking effects discussed above. To explain the salient 
points, it is easiest to initially consider a simpler field the- 
ory than the Q-matrix theory studied above, and adapt 
a classical line of reasoning from Ref. [p| to quantum field 
theories. Accordingly, we consider a scalar quantum field 
0(x, t) and an action 



S[(t)\^ / dx (^9^(A-H[0,V(/)]) 



(6) 



Here x = (x, r) comprises position x and imaginary time 
T, J dx = J dTi J dr, Ti. is a Hamiltonian density, and we 
use units such that h = k^ = 1. We will assume that 
S describes a phase transition from a disordered to an 
ordered phase, and will use a magnetic language, refer- 
ring to {(p) as 'magnetization'. Suppose that Ti. contains 
quenched disorder of random-mass type, and that we are 
in the nonmagnetic phase, (0) — 0. The key idea is to 
not integrate out the quenched disorder as a first step, as 
one does in a conventional treatment B, but rather to 
work with a particular disorder realization. Due to the 
quenched disorder there will be regions in space that en- 
ergetically favor local order, ((/)) 7^ 0, even though there is 
no global order. Deep inside the disordered phase these 



regions will be rare, but in an arbitrarily large system 
we will find arbitrarily large rare regions with a finite 
probability. The action S will then have static saddle- 
point solutions <E>(x) that have a nonvanishing value of 
the magnetization only in the rare regions. Let there be 
N such rare regions and associated local blobs of magne- 
tization or LMs. Then we can actually construct 2^ such 
saddle points, which differ only by the way the sign of the 
magnetization is distributed among the LMs. Since the 
LMs are far apart, the energy differences between these 
2^ saddle points will be small. In expanding about the 
saddle points, we therefore have no reason to prefer one of 
them over any of the others. Furthermore, since the LMs 
are self-generated by the system, albeit in response to the 
quenched potential, we assume that they are in thermal 
equilibrium with all other degrees of freedom as well as 
with each other. To calculate the partition function Z it 
is therefore necessary to take into account fluctuations in 
the vicinity of each of the 2^ saddle points M : 



Zk I D[ip\ e'^I'^1 I D[^ exp(^ f dx^'^{x) 

X exp ( / dx <I)(x) (fi'^ix) ] . (9) 



2~ 

Z^J2j ^[^] exp (-^ic&f'^) 



f] 



(7) 



Here J D[ip] denotes an integration over small fluctua- 
tions If in the vicinity of each of the saddle points. Notice 
that this restriction to small fluctuations is necessary in 
order to avoid double counting. Conversely, if we could 
perform the integral over the fluctuations exactly, then 
it would be sufficient to expand about one of the saddle 
points. In practice, however, one is restricted to a per- 
turbative evaluation of the functional integral, and Eq. 
(0) is a good approximation |ll[ |. 

We now consider the thermodynamic limit. Then the 
discrete set of 2^ saddle points turns into a saddle- 
point manifold Al($) that needs to be integrated over. 
Splitting off the saddle-point part of the action, 5'[(/)] — 
5[$] + AS[<^,ip], we have 

£>[$] P[$] /i)[(^] exp (-A5 [$,(/?]) , (8a) 

with the probability distribution P given by 

P[$] =5($) exp(-;^ /dx7^[$,V$] j . (8b) 

Here S denotes the support of the saddle-point manifold 
Ai. Notice the factor of l/T in the exponent, which 
results from the static nature of the saddle points |l2| . 

In general it is not possible to determine P[^] explic- 
itly. However, if we perform the $ integration by means 
of a cumulant expansion, the most relevant term in the 
effective action will be the one that results from the term 
quadratic in 7i[<I'] and the linear coupling between $ and 
ip'^ in AS. To obtain the most relevant term in the ef- 
fective theory for the fluctuations ip, we thus can write, 
with w > a number ]13] , 



Equation (ph is the partition function one would ob- 
tain by expanding perturbatively about just one of the 
saddle points, with static, annealed disorder appearing 
in addition to the quenched disorder still contained in 
S[ip\. The annealed disorder is governed by a Gaussian 
distribution whose variance is proportional to T. This 
property reflects the fact that the annealed disorder, as 
classical degrees of freedom in equilibrium with the rest 
of the system, must come with a Boltzmann weight, and 
it is the reason for the factors of T in Eqs. (||). 

Let us now explain how these arguments can be applied 
to the Q-field theory of interacting electrons to arrive at 
the action, Eqs. (|l|, g). The magnetization is propor- 
tional to the expectation value (tr (T3 (g) Si) Q(x)) 0, and 
in the presence of quenched disorder that favors the for- 
mation of magnetic LMs, the exact fermionic theory that 
underlies the NLctM |j] allows for saddle-point solutions 
where these components of Q are locally nonzero and 
play the role of the field $ above. This is in addition 
to a globally nonzero (tr (tq (8) sq) (5(x)) which reflects a 
nonvanishing DOS. By following the above reasoning for 
a scalar field, and going through the derivation of the 
sigma-model again, one obtains Eqs. (0,|^). 

We conclude with several remarks. First, we emphasize 
that we have studied a simplified model, neglecting both 
the Cooper channel and the electron-electron interaction. 
The latter point requires some clarification. In order to 
generate the annealed disorder from LMs, some interac- 
tion is necessary, (1) for local magnetic order to develop, 
and (2) in order for our canonical averaging over the sad- 
dle points to make physical sense. A truly noninteracting 
system would not sample all of these field configurations. 
Put differently, interactions make the energy barriers be- 
tween the saddle points, which are infinite in a noninter- 
acting system, finite and thus allow for an equilibration of 
the saddle-point degrees of freedom [p|Jll|] . We have sim- 
plified our model by assuming points (1) and (2) above to 
be the only effect of the interactions. Of course, if the an- 
nealed disorder were due to some other mechanism, then 
our results would also apply to strictly noninteracting 
electrons. Clearly, one can study generalizations of our 
model. In addition to adding an explicit interaction term, 
one can restore the Cooper channel, which will make the 
FP we found compete with the ordinary localization FP 
that also occurs at one-loop order. In systems with time 
reversal symmetry, one then expects the MIT studied 
here to get preempted by a localization transition if the 
bare dimcnsionless mass M2/ H is smaller than a num- 
ber of 0(1). It would also be interesting to consider the 
present model to 2-loop order to see whether the diffu- 
sion coefficient will still not be renormalized (apart from 



the 'difFuson' localization contributions that will appear 
at that order), and whether the line of FPs gives way to 
a more conventional FP structure. These questions will 
be considered in the future. 

Second, we point out that the strong effects of annealed 
disorder we found are characteristic of quantum statisti- 
cal mechanics. In a classical scalar field theory, the lead- 
ing term in the action generated upon integrating out an- 
nealed disorder is of the form (see Eq. (^)) — / dx </J^(x). 
It thus has the same form as the ordinary (/s -term and is 
in general not very interesting (although it can lead, e.g., 
to a first order phase transition). In a quantum system, 
on the other hand, integrating out the annealed disor- 
der yields — / dx / dr dr' ip'^ (x, r) ip'^ (x, r') , which has a 
different time structure than the usual ip^ term. It is 
the extra time integral that makes the annealed disorder 
term more relevant than in the classical case. 

Third, we come back to the fact that the variance of the 
Gaussian distribution for the annealed disorder is linear 
in T. If one used a Gaussian distribution with a temper- 
ature independent width, one would encounter factors of 
1/T in perturbation theory that force one to scale the an- 
nealed disorder strength with T to obtain a meaningful 
theory. Annealed disorder with an unbounded distribu- 
tion and a finite variance at T = is unphysical, since 
it allows the system to lower its energy arbitrarily far by 
digging itself a deeper and deeper trough. The necessity 
of the factor of T was realized in Ref. pfl] , but its origin 
was not recognized |l3]. 

Finally, let us explain why annealed disorder leads 
to a critical dn/dfi, while quenched disorder without 
electron-electron interactions does not. To see this, we 
realize that annealed disorder essentially means poten- 
tial troughs that are somewhat flexible, i.e. they adjust 
in response to the electrons. Let the sytem be in equi- 
librium at some value of the chemical potential fi, and 
change /i slightly. Then the flexible potential will adjust, 
and as a result fewer electrons will have to flow out of or 
into the grand canonical reservoir than would be the case 
in the absence of annealed disorder. This explains why 
there is a correction to dn/dfi in perturbation theory. 
Furthermore, the diffusive dynamics of the electrons lead 
to this correction being a frequency-momentum integral 
over diffusion propagators, which is logarithmically sin- 
gular in 2-d. In d = 2 + e this leads to a critical dn/dfj,, 
as it happens with other quantities that are singular in 
perturbation theory in 2-d. This is the only known mech- 
anism for a critical dn/dfi at a MIT in low-dimensional 
systems Wm . The recent observation of a critical dn/dfi 
at a 2-d MIT |jl3] is therefore very interesting in this con- 
text, even though our current theory does not describe a 
MIT in d = 2. 
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